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Simultaneous single-pulse observations of radio pulsars 
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Abstract. In this paper, we present a study of orthogonal polarization modes in the radio emission of PSR 
B1133+16, conducted within the frame of simultaneous, multi-frequency, single-pulse observations. Simultaneously 
observing at two frequencies (1.41 GHz and 4.85 GHz) provides the means to study the bandwidth of polarization 
features such as the polarization position angle. We find two main results. First, that there is a high degree of 
correlation between the polarization modes at the two frequencies. Secondly, the modes occur more equally and 
the fractional linear polarization decreases towards higher frequencies. We discuss this frequency evolution and 
propose propagation effects in the pulsar magnetosphere as its origin. 

Key words, pulsars: individual: PSR B 1133+ 16 - polarization 



1. Introduction 

Simultaneous multi-frequency observations of radio pul- 
sars in full polarization are proving to be a very pow- 
erful tool for investigating the p uzzling pulsar emission 
problem. In Karastergiou et al. ( 2001 , hereafter Paper 
I), we provided a detailed introduction into the history 
of multi-frequency observations and the ongoing efforts 
to construct an emission theory which will be both self- 
consistent and consistent with the abundance of obser- 
vational information from pulsars. We also laid out the 
guidelines for simultaneous polarimetric observations be- 
tween different telescopes by presenting our time-aligning 
and re-binning techniques. 

Observations of the Vela pulsar by Radhakrishnan & 



thus the linear polarization occurred in 'orthogonally po- 
larized modes' (OPM). The existence of OPM was a key 
discovery (Backer & Rankin 1980). By using simple argu- 
ments it can be shown that the presence of OPM acting 
together in the emission process can reduce the fraction of 
polarization observed (Stinebring et al. 1984, McKinnon 
1997, McKinnon & Stinebring 1998). Subsequently, single 
pulse polarization studies by Gil & Lyne (1995) demon- 
strated that OPM was responsible for the low percentage 
polarization and very complicated PA variations in the 
integrated profile of PSR B0329+54, in which two simple 
OPM modes existed. 



Melrose & Stoneham (1977) and Barnard & Arons 



Cooke (1969) showed the radiation to be highly linearly 
polarized and to be emitted from near the magnetic pole. 
The observed position angle (PA) reflects the instanta- 
neous projection of the dipolar magnetic field lines on the 
plane of the sky and thus changes with the rotation of the 
star (the so-called rotating vector model, RVM). Later, 
Manchester et al. (1975) demonstrated that the percent- 
age polarization could vary between and 100% across 
pulse longitude, that significant circular polarization was 
often present and, perhaps most surprisingly, that signif- 
icantly different PAs were seen at the same pulse longi- 
tude. These PA jumps were often very close to 90° and 



( 1986 ) proposed a model for the OPM in pulsars by con- 
sidering propagation effects (refraction) which cause rays 
in the two natural modes of the plasma to split. This has 
been further developed by McKinnon (1997) who showed 
that a model involving birefringence of the two modes can 
explain the observed pulse width and fractional polariza- 
tion changes with frequency. 

In this paper, we shall focus on the polarization be- 
haviour of PSR B1133+16 and use the results to draw 
some general conclusions about the pulsar emission pro- 
cesses. 
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Fig. 2. PA-PA plot for phase bin 50. The distribution of 
the PA values for each frequency is also plotted. 



Fig. 1. The total power profile of PSR B1133+16 is shown 
together with the linear polarization at two frequencies; 
the dotted line corresponds to 1.41 GHz and the dashed 
line to 4.85 GHz [TOP]. The grey-scale PA histograms 
follow, at 4.85 GHz [MIDDLE] and 1.41 GHz [BOTTOM]. 
The darker areas in the grey-scale histograms correspond 
to more frequent values of PA, as they occur on a single- 
pulse basis. The PA values are not absolute. 

2. Observations 

On 2000 January 4, we used the Effelsberg 100-m tele- 
scope at 4.85 GHz and the Lovell 76-m telescope in Jodrell 
Bank at 1.41 GHz to observe PSR B1133+16 (see von 



Hoensbroech et al. 1998, von Hoensbroech & Xilouris 1997 



for details on the Eff elsber g observing system and Paper 
I and Gould & Lyne 1998 for the Jodrell Bank system). 
Both telescopes were simultaneously on source for s=s 93 
minutes or 4778 pulse periods. For the purpose of this 
paper, we also make use of Jodre ll Bank data observed 
at 610 MHz (Gould & Lyne |199S| ), as well as Effelsberg 
2.69 GHz data (Paper I) from different epochs. The data 
from each telescope were calibrated and converted into 
the EPN- format (Lorimer et al. 1998] ) , after being aligned 
and re-binned to a common reference frame and temporal 
resolution (in this case 1.16 ms per bin), as described in 
Paper I. 

3. Results 

3.1. The behaviour of modes with frequency 

Fig. 1 shows the integrated pulse profile at 1.41 GHz along 
with the linear polarization at both frequencies. In this pa- 
per we ignore the circular polarization, which is generally 



low throughout the pulse. Fig. 1 also shows the PA distri- 
bution as a function of phase in a grey-scale representation 
at both 1.41 GHz and 4.85 GHz. It can be seen that at 
the peak of the first component two orthogonal modes are 
clearly present at both frequencies. In the second compo- 
nent one mode dominates at 1.4 GHz but two modes can 
be seen at the higher frequency. Integrated profiles like 
these have been used to question the broad-band nature 
of OPM (e.g Stinebring et al. 1984); for the first time, 
however, we can use the simultaneous nature of our data 
to determine this issue. 

At each of the two observed frequencies we compute 
the PA only if the flux density of the linear polarization 
exceeds twice the noise. We do this for each phase bin of 
every pulse; an example of one such bin near the peak of 
the first component is given in Fig ||. Under incoherent 
addition of the OPMs, the PA is only permitted to take 
one of two discrete values, exactly 90° apart. In this bin, 
it is clear that the orthogonal PA modes are present at 
both frequencies, although the spread in PA of each mode 
is somewhat larger than one would expect from a delta- 
function convolved with instrumental noise (Stinebring et 
al. 1984, McKinnon & Stinebring 1998,). Two clear re- 
sults are apparent in the figure. First, at 1.41 GHz the 
ratio of pulses in the two modes is 65:35 whereas at 4.85 
GHz this ratio is 54:46. The modes are occurring more 
equally at 4.85 GHz and the fractional linear polarization 
has also decreased. Secondly, all four possible combina- 
tions between the two modes at the two frequencies occur, 
i.e. four islands are visible in Fig. |^. Given the difference 
in the bimodal distributions at the two frequencies, such 
a distribution of points in the plot must arise, even in the 
case when there is no correlation between the PAs of the 
two frequencies. Hence, it is not readily clear from Fig. 2 
whether the cases of agreement and disagreement of OPM, 
represented by the four islands, suggest a high or low de- 
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Fig. 3. A sequence of phase bins from the leading component of the pulse. Each frame consists of a PA-PA histogram 
for the phase bin in the lower left corner, as identified in Fig. 1, and shows the deviation of the observed PA pairs 
from the null hypothesis, weighted by the calculated a (see explanation in text). 



gree of correlation between the frequencies. Obviously, the 
correlation is not perfect since the individual distributions 
are different — which is an important result already - but 
some degree of correlation may still exist. We have to test 
whether the distribution of points in the PA- PA plot can 
be attributed to some random combination of PA-values 
at different frequencies, or whether the knowledge of the 
polarization mode at one frequency allows one to predict 
the mode at the other to some extent. We test the degree 
of correlation by comparing the observed PA- PA distribu- 
tions to the null hypothesis that the two frequencies are 
completely uncorrelated. 

We re-bin our PA distributions into histograms of 12 
bins across the 180° range of PA values. By doing this 
we create a 2-dimensional 12 x 12 histogram Obs(i,j), 
which represents the number of occurrences of all possible 
PA-PA pairs as they are observed. Then we use Monte 
Carlo simulations (MCS) to randomly combine our two 
observed PA distributions. We repeat this process a large 
number of times, thus obtaining another 2-dimensional 
histogram Null(i,j), similar to Obs(i,j), where the popu- 
lation in each square bin is the mean result from all MCS 
iterations. At the same time, we can calculate the stan- 
dard deviation, a, representing the statistical fluctuations 
which we can expect in each square bin of the null hypoth- 
esis case. We now subtract the null hypothesis Null(i,j) 
from the observed 2-d histogram Obs(i,j), and compare 
this difference with the expected noise in each square bin, 
a. This will give us the significance in units of a that an 
observed number of points Obs(i,j) is not consistent with 
the null-hypothesis of no correlation. 

Such a null hypothesis test can only be applied reason- 
ably to bins where we observe bimodal PA distributions at 
both frequencies, as other cases cannot be distinguished 
from a random combination of PA values per se. Hence, 
in Fig. H we plot a 2-d PA- PA histogram only for the phase 
bins of the leading pulse component where the individual 
PA distributions are bimodal. The solid areas denote a 
difference between Obs(i,j) and Null(i,j) of more than 
3cr surrounded by a lighter coloured area marking differ- 
ences exceeding 2a. Hatched areas denote differences of 
less than — 3(7 surrounded again by lighter coloured area 
of —2(T significance. We can clearly see that when one 
mode occurs at one frequency, it also prefers to occur at 
the other frequency: the same modes combine much more 
times than the null hypothesis of no-correlation predicts 
(black areas on the x=y diagonal). 



In cases where the modes disagree, one mode does not 
necessarily favour a value at the other frequency that is 
different by exactly 90°. On the contrary, for most of the 
time the PA value at the other frequency seems to be 
completely unrelated, resulting in a distribution of points 
outside the diagonal islands that is mostly consistent with 
the null hypothesis. This is an interesting result with pos- 
sible implications on the widths of PA distributions at 
given frequencies. Only in few cases is a lack of points 
observed when compared to the null hypothesis (hatched 
areas), supporting an increased degree of correlation be- 
tween the frequencies. We note that Fig. || does not change 
qualitatively by using grids coarser than 12 x 12. 

The results in this section therefore show that there is 
a high degree of correlation between the two frequencies. 
The same mode tends to dominate at both frequencies and 
to this extent OPM is broadband in nature. At the same 
time, the modes occur more equally and the fractional 
polarization is lower at the higher frequency. This will be 
discussed in more detail below. 



3.2. The x and s parameters 



According to Stinebring et al. (|1984[ ), the emission we re- 
ceive at any pulse phase is the instantaneous, incoherent 
superposition of the two OPMs. The PA value is that of 
the strongest mode. Two important parameters can be 
constructed to investigate the strength and behaviour of 
the modes. 

The first parameter, x, represents the relative popu- 
larity of the PA of one mode with respect to the other. 
By studying the individual PA distributions for each fre- 
quency, it is possible to identify the OPMs across the 
phase bins of the pulse. We can identify the PA of each 
mode in each bin i by finding the mean PA (or tracing 
the PA swing when the mode seems absent) and identi- 
fying a PA range of « 30° each side of the mean, which 
corresponds to the typical width of the PA histograms, as 
shown in Fig. ||. Mi is the number of pulses where the PA 
corresponds to mode 1 and M2 the same for mode 2. For 
each of the frequencies we studied the occurrence of Mi 
and Mi at every pulse phase bin. 

From the values of M\ and M2 we have calculated for 
each phase bin, we construct x — ffSrffj- The value of x 
depends on the relative values of M\ and M-z- If they are 
equal, x will obviously be zero. On the other hand, if one 
of them is larger than the other, x will diverge from zero, 
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Fig. 4. The parameter x [solid line] and \s\ [dashed line] 
across the phase bins of PSR B1133+16 for four different 
frequencies. There is a general trend for both x and s to 
become closer to zero toward higher frequencies. 



1 and —1 being the extreme cases, where only one of the 
modes is present. 

The second parameter is s, which represents the rel- 
ative strengths of the modes and varies from —1 to 1 in 
a similar fashion to x. s is defined as s = i 1 ~ s a 2 , where 



each phase bin is assumed to be the incoherent sum of 
two 100% orthogonally polarized modes with individual 



flux densities Si and S2 (Stinebring et al. 1984). In the 
case where the modes have a negligible fraction of circular 
polarization (as is the case here), then |s| = L/I, where L 
is the linear polarization and / the total power. 

We used the simultaneously observed data to study 
the behaviour of x and s with frequency. We noticed that, 
although they are both phase dependent, they generally 
seem to have values closer to zero at 4.85 GHz than at 
1.41 GHz (Fig. |). We then used other observations of 
PSR B1133+16 to see if this behaviour is consistent at 
other frequencies. The Jodrcll Bank 610 MHz data and the 
Effclsberg 2.69 GHz data, which are of comparable length 
to the simultaneous data, clearly show this remarkable 
trend. In other words, the times M2 or Mi determine the 
PA become more equal towards higher frequencies, while 
at the same time mode 2 is becoming closer in strength to 
mode 1 (s closer to 0). 



4. Discussion 

In the case of PSR B1133+16 we have found the following 
trends of the OPMs: 

— OPMs have a spectral dependence: the polarized in- 
tensities of the modes (s) and their frequency of oc- 
currence (x) are becoming comparable at higher fre- 
quencies; 

— we observe phase bins that simultaneously show a high 
value of x and a low value of s; 



— despite the spectral evolution, we find a strong con- 
nection between the two frequencies: the same polar- 
ization modes prefer to appear together; 

— different modes, however, do appear at the two fre- 
quencies, implying a degree of de-correlation. 

The first two trends are also supported by the 610 MHz 
and 2.69 GHz data. 

How well do these observational trends fit in with the 
model of superposed orthogonal modes? A high value of 
x implies that a given mode dominates most of the time. 
A low value of s implies that the difference in the average 
strengths of these two modes are rather small. These two 
facts can only be reconciled by postulating that the ampli- 
tudes of the two modes are highly correlated (McKinnon & 
Stinebring 1998). A physical mechanism which can induce 
mode splitting with a high degree of amplitude correlation 
is birefringence and this means of mode creation in the rel- 
ativistic pulsar plasma was proposed by Barnard & Arons 
(1986). 

We also observe that, at least in this pulsar, x decreases 
more rapidly with frequency than does s. This clearly im- 
plies that the correlation coefficient is decreasing with in- 
creasing frequency even though the relative strengths of 
the modes change rather little. One observational conse- 
quence of this is that we would expect that the distribution 
of linear polarization in the single pulses has a rather low 
dispersion when the correlation coefficient is high, com- 
pared to when the correlation coefficient is low. This is 
indeed seen in the data; the distribution of linear polar- 
ization is broader in the single pulses observed at 4.85 
GHz. 

We therefore surmise that the relative strength of the 
modes is set by the polarization of the input beam to the 
bi-refringent medium and that this is largely constant with 
frequency. The frequency trend in the observed value of 
s arises from the refractive index of the modes being fre- 
quency dependent - at some 'critical' frequency and above, 
the beams overlap, causing a decrease in s. In this pulsar, 
s changes rather little with frequency, as do the mode- 
separated pulse widths between 1.41 GHz and 4.85 GHz 
and presumably the critical frequency lies near 1.4 GHz. In 
most other pulsars, however, the critical frequency seems 
to lie above 5 GHz (Xilouris et al. 1996 ). 

The third and fourth conclusions we reach, together 
with the fact that total power is well correlated in our 
data (as was the case for PSR B0329+54 in Paper I), 
also point towards propagation effects. They support the 
idea that OPMs, which could be created by the splitting 
of the two natural orthogonal polarization modes in the 
magnetosphere by a propagation effect like birefringence, 
endure further propagation effects that change their rel- 
ative strengths. One explanation could be that the ab- 
sorption (e.g. cyclotron absorption) of one or other of the 
modes is frequency dependent with the higher frequencies 
more affected than lower frequencies. A change in relative 
strength could switch the dominant mode. The obvious 
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manifestations of this are the clearly seen OPM disagree- 
ments at the two frequencies. 

A further observational test of propagation effects and 
the superposed mode ideas comes from the circular polar- 
ization. One expects that a given mode is associated with 
a given sign of circular polarization and that there is there- 
fore a strong correlation between the dominant mode and 
the sign of circular polarization in the single pulse data. 
Given that we see a high correlation between mode occur- 
rences at different frequencies, will we also see a correla- 
tion between the handedness of the circular polarization? 
Simultaneous, multi- frequency, single pulse observations 
are ideal as a test to identify propagation effects at work 
on the radiation traveling through the plasma in the pul- 
sar magnctosphere. 
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